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Abstract 

First-principles  calculations  within  the  spin-polarized  generalized  gradient  approximation  framework  are  performed  on  rhombohedral  and  its 
Jahn-Teller  distorted  monoclinic  LiMnC>2  with  different  spin  configurations.  It  is  found  that  the  Jahn-Teller  (JT)  effect  is  driven  by  high- spin  Mn3+ 
ion  and  it  should  disappear  in  low- spin  state.  With  JT  distortion,  the  initially  degenerate  eg  states  is  split  by  a  gap  of  687  meV,  which  responds  for 
the  semiconduction  of  monoclinic  LiMn02-  Based  on  the  analyses  of  the  changes  induced  by  JT  distortion  in  the  crystal,  electronic  structures  and 
chemical  bondings,  approaches  to  suppress  the  JT  effect  and  synthesize  rhombohedral  LiMn02  are  suggested.  At  last,  the  JT  effect  is  decomposed 
into  an  electronic  and  an  elastic  term. 
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1.  Introduction 

Because  monoclinic  LiMnC>2  (space  group  C2/ra,  here¬ 
after  denoted  as  ra-LiMn02)  has  a  high  theoretical  capacity 
(285mAhg_1)  and  Mn  is  abundant,  inexpensive  and  environ¬ 
mentally  benign,  intense  expectation  has  been  put  on  this  mate¬ 
rial  to  replace  LiCoC>2,  the  current  standard  cathode  material 
in  commercial  Li-ion  rechargeable  batteries  since  it  was  first 
obtained  by  ion  exchange  from  a-NaMnC>2  in  1996  [1,2].  This 
compound  has  a  layered  structure  where  Li+,  Mn3+  and  02~ 
locate  in  their  own  triangular  planes,  respectively,  and  these 
planes  stack  with  a  sequence  as  Li-O-Mn-O-Li.  The  O2-  sub¬ 
lattice  is  fee  and  simply  provides  a  framework  with  Li+  and  Mn3+ 
filling  all  the  interstices  of  its  edge-shared  octahedrons.  Such  a 
layered  structure  is  similar  to  that  of  LiCoC>2  and  has  superior  Li 
de-intercalation  cycling  properties  than  that  of  its  stable  form, 
orthorhombic  LiMn02  ( Pmnm ,  6>-LiMn02)  [1-5].  Experimen¬ 
tally,  a  large  initial  capacity  of  more  than  270  mAh  g~ 1  has  been 
measured  [1]. 


*  Corresponding  author.  Tel.:  +86  431  5168446;  fax:  +86  431  5168446. 
E-mail  address:  gchen@jlu.edu.cn  (G.  Chen). 

0378-7753/$  -  see  front  matter  ©  2005  Elsevier  B.V.  All  rights  reserved, 
doi :  10. 1016/j  .jpowsour.2005 . 1 0.044 


However,  m-LiMn02  suffers  from  severe  capacity  fading, 
which  is  surely  ascribable  to  the  irreversible  transition  to  a 
spinel-like  structure  and  is  the  primary  obstacle  to  its  appli¬ 
cations  in  practice  [3-5].  Such  a  structural  instability  is  caused 
by  the  wide  and  discontinuous  variation  in  the  magnitude  of 
the  Jahn-Teller  (JT)  distortion  because  the  Mn  ions  are  cycled 
between  JT  active  high-spin  Mn3+  and  inactive  Mn4+  states  dur¬ 
ing  the  electrochemical  charge-discharge  process  [5,7].  There¬ 
fore,  JT  effect  plays  an  important  role  in  the  applications  of 
this  electrochemical  material  and  deserves  detailed  investiga¬ 
tion.  JT  effect  operates  due  to  partial-filled  electronic  sub¬ 
shell,  distorting  the  crystal  structure  to  eliminate  the  degen¬ 
eracy  until  the  increase  in  elastic  energy  halts  the  distortion 
[5-7]. 

Because  m-LiMnC>2  is  the  JT  distorted  form  of  rhombohe¬ 
dral  LiMnC>2  ( R3m ,  r-LiMn02),  we  carry  out  first-principles 
calculations  on  the  crystal  and  electronic  structures  of  both  m- 
and  r-LiMnC>2  with  different  spin  configurations  to  investigate 
the  JT  effect  in  ra-LiMn02.  Previous  reports  have  concluded 
that  generalized  gradient  approximation  (GGA)  achieves  greater 
accuracy  for  Mn-oxides  compared  with  local  density  approxima¬ 
tion  (LDA)  [8].  In  addition,  Mn3+  ions  in  both  the  structures  are 
magnetic.  Therefore,  spin-polarized  GGA-PW91  (GGS-PW91) 
[9]  is  employed  for  the  exchange-correlation  potential. 
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2.  Computational  details 

All  the  calculations  are  performed  with  the  CASTEP  soft¬ 
ware  [9].  The  initial  crystal  model  of  m-LiMn02  is  built  using 
the  experimental  lattice  parameters  and  ionic  positions  [1]  but 
neglecting  the  ~10%  Li/Mn  site  disorder,  while  hypothetical 
r-LiMn02  has  isostructure  to  that  of  LiCo02  [10]  but  sub¬ 
stituting  all  Co  with  Mn.  Primitive  cells  containing  one  for¬ 
mula  unit  are  investigated.  In  an  octahedral  complex,  the  five 
d  orbitals  of  a  transition  metal  ion  will  be  split  into  a  set 
of  lower  three-fold  degenerate  t2g  states  and  upper  two-fold 
degenerate  eg  states  according  to  ligand  field  theory  [6,7]. 
Therefore,  Mn3+  (d4)  ion  can  have  three  spin  configurations 
as  high-spin  (t3CTeg),  low-spin  (t4CTeg  and  occupying  all  three  t2g 

levels)  or  non- spin  (t4CTeg  and  occupying  two  of  the  three  t2g 
levels)  states,  which  is  stimulated  by  fixing  the  spin,  unpaired 
electrons  in  the  cell,  as  4,  2  or  0  under  neutralizing  back¬ 
ground,  respectively.  Ultrasoft  pseudo-potentials  in  the  recip¬ 
rocal  space  are  employed  with  a  plane-wave  energy  cutoff  of 
400  eV.  Requested  k-point  spacing  is  set  to  0.05  A-1,  which 
corresponds  to  72  k-points  for  m-LiMn02  and  60  k-points 
for  r-LiMn02  in  the  irreducible  Brillouin  zone  generated  by 
the  Monkhorst-Pack  scheme.  Geometry  optimization  is  first 
performed  to  fully  relax  the  structures  until  self-consistent 
field  (SCF)  convergence  per  atom,  tolerances  for  total  energy, 
root-mean- square  (RMS)  displacement  of  atoms,  RMS  force 
on  atoms  and  RMS  stress  tensor  are  less  than  2  x  10“6  eV, 
2xlO-5eV,  0.001  A,  0.05  eV/A  and  0.1  GPa,  respectively. 
Then,  the  relaxed  structures  were  used  to  calculate  the  total 
energies  and  the  electronic  structures  with  the  same  energy  cut¬ 
off  and  k-points  as  those  the  geometrical  optimizations  used. 
Because  m-LiMn02  shows  a  spin-glass  ground  state  where 
short-range  anti-ferromagnetic  (AF)  and  ferromagnetic  (FM) 
interactions  are  coexistent  randomly  [3],  for  the  restriction  of 
invoked  periodic-boundary  conditions  during  the  calculations, 
only  FM  ordering  was  considered  with  disregarding  the  AF 
cases. 


3.  Results  and  discussion 

The  calculated  relative  total  energies  (A F),  magnetic 
moments  at  each  Mn  site  (M)  and  lattice  parameters  (a,  b,  c 
and  P)  of  the  relaxed  crystal  structures  together  with  the  exper¬ 
imental  and  theoretical  data  are  listed  in  Table  1.  As  proven 
by  M  obviously,  the  high-spin,  low-spin  and  non-spin  states 
have  been  reproduced  successfully  (the  reduction  of  M  from 
the  nominal  Hund’s  rule  value  for  high- spin  state,  i.e.  M  are  not 
equal  to  but  less  than  4.0,  will  be  discussed  below  in  the  sec¬ 
tion  of  partial  density  of  states).  From  the  viewpoint  of  energy, 
non- spin  state  is  instable  distinctly  to  low- spin  and  high- spin 
states  for  both  symmetries,  so  this  spin  configuration  will  not 
be  considered  hereafter.  For  high-spin  state,  because  of  the 
energy  benefit  as  large  as  175  me V,  the  JT  effect  will  usually 
take  place  during  the  progress  of  synthesis  where  the  experi¬ 
mental  condition  can  supply  the  additional  energy  to  overcome 
the  energy  barrier  to  the  transition,  and  thereby  m-FiMn02  but 
not  r-FiMnC>2  is  obtained.  As  to  low-spin  state,  situation  is  dis¬ 
cussible  because  the  energy  benefit  from  JT  distortion  is  too 
small  (only  8meV).  Such  small  value  may  locate  in  the  calcu¬ 
lation  error  and  more  importantly,  according  to  the  argument 
of  Ceder  et  al.,  a  small  enough  (less  than  30meV)  formation 
enthalpy  allows  for  entropy- driven  mixing  [  1 1  ] .  In  the  same  way, 
low-spin  state  of  r-FiMnC>2  can  be  stabilized  to  high-spin  state 
with  an  energy  barrier  of  only  1 1  meV.  Experimentally,  high- 
spin  state  for  m-FiMn02  and  low-spin  state  for  r-FiMnC>2  also 
reproduce  compatible  lattice  parameters  with  measured  data  (see 
Table  1)  [1,12].  Because  pure  r-FiMnC>2  is  still  hypothetical,  the 
Cr-doped  experimental  data  are  used  here  [12].  Theoretically, 
other  groups  have  found  a  low-spin  Mn3+  ion  in  r-FiMnC>2,  too 
[5,7,8].  Therefore,  it  is  concluded  that  the  JT  distortion  in  lay¬ 
ered  FiMn02  is  driven  by  the  high-spin  Mn3+  ion,  resulting  in 
m-FiMn02.  While  for  low-spin  state,  JT  distortion  is  absent  and 
r-FiMnC>2  should  be  obtained. 

Table  1  also  gives  the  nearest  Mn— Mn  distances,  which  are 
equal  to  the  lattice  parameter  b,  cell  volumes  (V)  per  formula 


Table  1 

Calculated  magnetic  moments  at  each  Mn  sites  (M),  relative  total  energies  (A E)  and  cell  volumes  (V)  per  formula  unit,  optimized  lattice  parameters  (a,  b,  c  and  ft), 
Mn— O  bond  lengths  (Lmh-o)  and  Jahn-Teller  distortion  degree  (JT)  of  r-  and  m-LiMn02  under  different  spin  configurations  (Sfix)  together  with  experimental  and 
theoretical  data 


SflX 

Cell 

Mmu  (p^b) 

A  E  (meV) 

a  (A)  b  (A) 

c  (A) 

PO 

V(A3) 

o 

TMn-0  (A) 

JT  (%) 

4 

r-LiMn02 

3.80 

0 

a  =  b  =  3.0045 

14.7392 

90 

38.4085 

2.0665  x  6 

0 

m-LiMn02 

3.72 

-175 

5.5497  2.8623 

5.4365 

116.575 

38.6167 

2.3707  x  2, 
1.9679  x4 

20.02 

2 

r-LiMn02 

2.00 

11 

a-b-  2.8664 

14.6002 

90 

34.6296 

2.0021  x  6 

0 

m-LiMn02 

2.00 

3 

4.9818  2.8628 

5.1237 

109.315 

34.4801 

2.0023  x  2, 
2.0014  x  4 

0.04 

0 

r-LiMnC>2 

0.00 

630 

a  =  b  =  2M09 

14.4659 

90 

33.7034 

1.9916x6 

0 

m-LiMn02 

0.00 

659 

5.0030  2.8675 

5.0760 

109.684 

34.2826 

2.0007  x  2, 
1.9896  x4 

0.56 

Experimental 

r-LiMn02  [12] 

— 

— 

a  =  b  =  2.8799 

14.3799 

90 

34.4287 

1.9718x6 

0 

m-LiMn02  [1] 

— 

— 

5.4387  2.8086 

5.3878 

116.006 

36.9837 

2.3079  x  2, 
1.9211  x  4 

18.87 

Theoretical 

r-LiMn02  [8] 

~2.0 

— 

a  —  b  —  2.82 

14.27 

90 

32.76 

1.97x6 

0 

m-LiMn02  [8] 

~4.0 

" 

5.54  2.82 

5.44 

116 

38.19 

2.34  x  2, 

1.92  x  4 

20.39 
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unit  and  Mn— O  bond  lengths  (LMn-o  )•  It  is  noted  that  with 
the  JT  distortion,  LMn-o  in  the  regular  MnC>6  octahedron  have 
split  into  two  sets:  two  are  elongated  but  the  other  four  short¬ 
ened.  Defining  JT  distortion  degree  as  the  difference  between 
the  longest  and  shortest  Mn— O  distances  divided  by  the  average 
Mn— O  distance  [3,8],  large  JT  distortion  (20.02%)  are  detected 
for  high-spin  state,  while  for  low-spin  state,  it  is  too  small  (only 
0.04%)  to  be  detected  experimentally.  Furthermore,  the  JT  tran¬ 
sition  temperature  ( Tt )  and  the  JT  distortion  degree  da  (here  c 
and  a  refer  to  the  inner  parameters  of  MnC>6  octahedron  but  not 
the  lattice  parameters)  have  a  moderate  correlation  to  the  elec¬ 
tronic  structure  of  the  MnC>6  octahedron  because  of  the  thermal 
excitation  of  electrons  over  the  energy  gap  A  between  the  &zi 
and  dr2_y2  orbital  of  eg  states  in  the  JT  distorted,  tetragonal 
symmetry  [13], 

A  ~  kg  T{  oc  (c/a  —  1)  (1) 

For  c!a- 1.15  and  Lt=1400K  [13],  we  can  extrapolate 
Tt  =  1910  K  for  high-spin  state  with  da  being  1.20.  Obviously, 
m-LiMn02  is  molten  at  this  temperature,  which  implies  that 
in  high-spin  state  the  JT  distortion  always  exists  in  the  crys¬ 
talline  structure  and  non-JT  distorted  r-LiMn02  cannot  be 
synthesized,  agreeing  well  with  the  experiments  [1-7].  For 
low-spin  state,  Lt  =  4K,  which  means  that  JT  distortion  will 
disappear  at  the  temperature  higher  than  4K  and  the  mate¬ 
rial  will  adopt  the  rhombohedral  structure.  That  is  to  say, 
in  order  to  obtain  r-LiMn02,  one  should  manage  to  keep 
the  material  in  the  low-spin  state,  reinforcing  the  conclusion 
above. 

In  order  to  compare  the  lattice  parameters  in  the  same  scale, 
monoclinic  C2/m  symmetry  can  be  indexed  by  the  equivalent 
lattice  parameters  am  and  cm  like  as  in  the  R3m  structure:  am  is 
calculated  from  the  area  of  the  triangular  net  (S)  and  cm  from  the 
interlayer  distance  (D)  based  on  the  relation  between  these  two 
structures  [4].  In  both  the  structures,  the  Li+,  Mn3+  and  O2-  ions 
form  respective  triangular  layers,  and  these  layers  stack  with 
a  sequence  as  Mn-O-Li-O-Mn.  Differently,  for  R3m  struc¬ 
ture,  the  triangular  net  of  the  layer  is  equilateral  with  the  side 
indexed  as  the  a  parameter,  therefore,  S  =  1/2  x  a  x  V3a/2. 
Because  its  unit  cell  contains  three  periods  of  the  sequence 
which  is  stacked  vertically  and  the  total  interlayer  distance  is 
indexed  by  the  c  parameter,  D  =  1/3  x  c.  While  for  C2/m  struc¬ 
ture,  JT  distortion  has  deformed  the  equilateral  triangular  net 
of  the  layers  to  be  isosceles  triangle  and  the  unit  cell  contains 
only  one  period  of  the  sequence  which  is  not  stacked  verti¬ 
cally  any  longer,  making  S  =  1/2  x  b  x  a/2  and  D  =  c  x  sin  f$. 
Therefore, 

/  a  x  b 

am  =  VZT 

Values  of  am  =  2.9849  A,  cm=  14.5864  A  and  am  =  2.8695  A, 

o 

cm  =  14.5059  A  are  obtained  for  high-spin  and  low-spin  states 
of  m-LiMn02,  respectively,  where  am  is  larger  and  cm  is  smaller 
than  those  for  r-LiMn02  with  the  same  spin  configuration  (see 
Table  1 ).  It  is  interesting  to  note  that  the  relaxation  is  not  achieved 
by  simply  dispersing  all  the  ions  in  the  layers  but  by  making  two 


of  the  six  coordinates  around  one  ions  closer  along  the  b  axis 
and  the  other  four  much  farther,  just  like  the  change  in  the  MnC>6 
octahedron  as  revealed  by  the  nearest  Mn— Mn  distances  b  (see 
Table  1).  At  the  same  time,  V  increases  in  high- spin  state  but 
decreases  in  low- spin  states.  These  changes  with  JT  distortion 
imply  that  JT  distortion  has  compressed  the  structure  perpen¬ 
dicularly  to  the  triangular  layers  but  relaxed  the  structure  in  the 
layers,  and  the  compression  dominates  in  low-spin  state  while 
the  relaxation  dominates  in  high-spin  state.  Therefore,  in  order  to 
suppress  the  JT  effect,  endeavors  should  be  introduced  to  weaken 
the  interaction  between  the  layers  but  strengthen  the  interaction 
in  the  layers.  Such  a  finding  can  answer  for  the  experimental 
fact  that  Co-doped  r-LiMn02  can  be  synthesized  and  it  tends  to 
form  non- stoichiometric  materials  [5],  because  the  absence  of 
Li  can  enhance  the  electrostatic  repulsion  between  the  nearest 
0—0  layers. 

As  described  in  the  introduction,  Mishra  and  Ceder  [8]  have 
performed  similar  theoretical  calculations  on  the  material  with 
special  emphasis  on  the  ground  states.  Their  results  for  FM  solu¬ 
tion  with  GGS  method  [8]  are  quoted  in  the  last  two  rows  of 
Table  1.  They  reported  that  the  low-spin  and  high-spin  states 
are  the  ground  state  of  r-  and  m-LiMn02,  respectively,  con¬ 
sistent  with  our  conclusion.  Nevertheless,  it  is  noticeable  that 
under  the  same  spin  states  our  values  for  the  lattice  param¬ 
eters,  V  and  LMn-o  are  all  larger  than  theirs,  which  is  more 
distinct  for  r-LiMn02.  We  ascribe  it  to  the  full  relaxation 
of  our  structures  because  of  the  severe  computational  condi¬ 
tions  used  during  the  geometry  optimization.  Obviously,  full 
relaxation  is  very  important  and  necessary  for  hypothetical 
r-LiMn02. 

Figs.  1  and  2  show  the  band  structures  around  the  Fermi  level 
(Lp)  for  high-spin  and  low-spin  states  before  and  after  JT  dis¬ 
tortion,  respectively.  Lp  is  set  as  OeV,  and  the  spin-up/down 
sub-bands  are  plotted  separately  as  denoted  by  “up”/“down” 
in  the  figures.  Fig.  3  is  the  corresponding  density  of  states 
(DOS).  According  to  the  schematic  picture  for  the  band  struc¬ 
tures  of  octahedrally  coordinated  solids  [14],  the  top  five  bands 
are  assigned  to  Mn-3d  bands  and  the  bottom  six  are  0-2p 
bands,  because  they  mainly  consist  of  Mn-3d  and  0-2p  states, 
respectively.  An  unoccupied  bands  located  around  lOeV  is 
mainly  made  up  of  Li-2s  and  Mn-4s,  4p  states,  and  the  valence 
bands  at  about  —18  and  — ■ 45  eV  are  0-2s  and  Li- Is  states, 
respectively  (not  shown  in  the  figures).  Such  assignments  are 
proven  directly  by  the  partial  DOS  as  shown  in  Fig.  3.  Because 
of  the  ligand  field  effect,  the  Mn-3d  bands  have  split  into 
upper  two  eg  bands  and  lower  three  t2g  bands  by  a  direct 
gap  (zA0). 

For  high-spin  state,  the  ground  state  of  m-LiMn02  as  proven 
above,  one  can  note  that  the  most  distinct  changes  in  the  bands 
structures  with  JT  distortion  are  that  the  two  eg  bands  have  been 
split  by  a  gap  of  687  me  V  which  responds  for  the  semiconduction 
of  this  material,  and  Ao  has  decreased.  Such  a  splitting  elimi¬ 
nates  the  degeneracy  in  half-filled  eg  bands  and  lowers  the  energy 
of  the  electrons,  thereby  the  system’s  energy.  The  later  cannot 
be  read  out  directly  from  the  bands  because  the  Lp  is  fixed  at 
0  eV.  For  comparison  these  bands  on  an  “absolute”  energy  scale, 
we  can  align  the  Li- Is  (or  0-2s)  states  based  on  the  assumption 


1/2 


cm  =  3c  x  sin 


(2) 
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Fig.  1.  Band  structures  for  high-spin  state  before  (a)  and  after  (b)  JT  distortion. 


that  these  states  are  core  states  and  should  not  be  affected  by  the 
different  environments  and  therefore  can  be  used  as  a  reference 
[15].  In  order  to  operate  this  alignment,  it  is  found  that  the  E p 
for  m-LiMn02  should  shift  down  by  1.041  eV  (see  Fig.  3).  That 
is  to  say,  eg  bands  have  shifted  down  about  1.041  eV  with  JT 
distortion.  However,  the  level  of  t2g  bands  remains  almost  fixed 
because  of  its  non-bonding  character.  Therefore,  it  is  inevitable 
that  Aq  decreases  dramatically.  The  decrease  in  Ao  also  con¬ 
sists  to  the  stretching  of  the  MnC>6  octahedron  caused  by  JT 
distortion,  because  Ao  characterizes  the  strength  of  the  ligand 
field  and  compressing  the  Mn— O  bond  length  would  lead  to  its 


(a)  r-LiMn02 


increase  [14,16].  When  the  compressing  is  considerable  enough 
to  split  these  two  bands  to  the  point  where  the  spin-up  eg  levels 
rise  above  the  lowest  spin-down  t2g  level,  low-spin  state  would 
be  obtained  from  the  high- spin  state  [16]. 

Actually,  it  is  just  the  splitting  of  eg  bands  that  reflects  the  JT 
effect  on  the  electronic  structures  and  characterizes  the  JT  dis¬ 
tortion  degree  [6,13,17].  Because  of  the  negligible  JT  distortion 
for  low-spin  state  as  revealed  above,  it  is  not  strange  that  the 
change  in  its  two  eg  bands  is  small  and  they  are  still  degenerate 
(see  Figs.  2  and  3).  Here,  it  is  necessary  to  point  out  that  JT 
effect  splits  the  eg  bands  with  the  same  spin  orientations  in  both 


Fig.  2.  Band  structures  for  low-spin  state  before  (a)  and  after  (b)  JT  distortion. 
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Energy  (eV) 


Fig.  3.  DOS  for  high-spin  state  before  (a)  and  after  (b)  JT  distortion  as  well  as  for  low-spin  state  before  (c)  and  after  (d)  JT  distortion.  Spin-up/down  states  are  plotted 
along  the  positive/negative  ordinate.  In  this  energy  range,  the  contributions  of  Li-s,  Mn-s,  Mn-p  and  O-s  are  small  and  are  not  shown.  The  energy  axes  have  been 
aligned  as  described  in  the  text. 


spin-up  and  spin-down  sub-bands,  while  the  splitting  between 
these  two  spin  species  is  caused  by  intra-ionic  exchange  [16,18] 
which  has  been  treated  wrongly  as  JT  effect  splitting  in  Ref. 
[7].  As  expected,  the  exchange  splitting  Z\exch  is  similar  before 
and  after  JT  distortion  because  of  the  similar  M  (see  Table  1). 
Taking  Aexch  as  the  energy  difference  between  the  midpoints 
of  the  spin-up  and  spin-down  t2g  bands  [18],  an  order  of  3.4 
and  1.8  eV  can  be  calculated  for  high-spin  and  low-spin  states, 
respectively. 

In  Fig.  3,  significant  mixing  of  spin-up  Mn-3d  states  in 
0-2p  bands  as  well  as  0-2p  states  in  eg  bands  can  be  noted 
in  all  compounds,  which  reveals  strong  hybridization  between 
Mn-3d  and  0-2p  states  and  the  bonding,  anti-bonding  and  non¬ 
bonding  characters  of  0-2p,  eg  and  t2g  bands,  respectively. 
This  hybridization  leads  to  three  consequences:  the  first  is  the 
exchange  splitting  of  0-2p  bands  because  of  the  exchange  split¬ 
ting  of  eg  states;  the  second  is  the  reduction  of  M  from  the 
nominal  Hund’s  rule  value  for  high- spin  state  (see  Table  1)  with 
the  residual  magnetic  moment  on  the  O  ions,  because  the  par¬ 
tially  filled  eg  states  are  spin  polarization;  the  third  is  the  partial 
ionization  of  O  ions  because  the  eg  states  are  not  fully  occu¬ 
pied.  Similar  results  are  reported  in  Singh’s  work  where  M  is 
only  3.45  jjib  and  Aexch  of  0-2p  bands  is  about  0.6  eV  [19]. 
At  the  same  time,  Li  ion  is  fully  ionized  as  revealed  by  the 
absence  of  Li-2s  states  near  E p.  In  order  to  keep  the  electric 
neutrality,  therefore,  Mn  ion  is  partially  ionized,  too.  That  is 
to  say,  the  number  of  3d  electrons  for  Mn  ion  is  larger  than 
4,  which  responds  to  the  deviation  from  the  normal  ionic  sce¬ 
nario  in  the  band  structures,  where  more  than  four  t2g  bands  lie 
below  for  low-spin  state  and  the  E for  high-spin  ra-LiMn02 
locates  above  the  spin-up  eg  bands  rather  than  in  the  middle  of 
them. 

The  large  degree  of  hybridization  also  indicates  strong  cova¬ 
lent  interaction  between  Mn  and  O  ions.  For  more  detailed  and 


deeply  discussion,  Mulliken  population  analysis  is  performed 
and  the  results  are  listed  in  Table  2.  The  strength  of  the  cova¬ 
lent/ionic  interaction  between  two  ions  is  characterized  by  a 
positive/negative  value  of  the  bond  overlap  populations  (OP). 
The  large  deviation  from  the  formal  charges  of  Mn  and  O  ions 
confirms  that  the  Mn— O  bonds  are  strongly  covalent,  which  is 
reinforced  by  the  large  positive  values  of  corresponding  OP. 
Comparably,  it  is  clear  the  Li  ions  has  been  absolutely  ionized 
and  the  Li— O  bonds  are  ionic  property.  Similar  pictures  occur 
in  other  lithium  transition  metal  oxides  and  it  is  found  vital 
to  the  materials’  electrochemical  performance  and  optimization 
[11,14,19-21].  For  high-spin  state,  it  is  found  that  the  JT  dis¬ 
tortion  has  increased  slightly  the  ionicities  of  Li  from  +1.03 
to  +1.04  but  decreased  the  ionicities  of  Mn  and  O  ions  from 
+0.61,  —0.82  to  +0.58,  —0.81,  respectively.  As  revealed  by  the 
increase  in  OP  values  and  decrease  in  the  bond  lengths,  the 
ionic  interaction  in  Li— O  bonds  and  the  covalent  interaction  in 
Mn— O  bonds  are  both  enhanced,  which  implies  that  the  Li  ion 
also  participates  in  the  JT  distortion  [22]  and  the  average  inter¬ 
calation  voltage  of  m-LiMn02  will  increase  [23].  Additionally, 
it  reveals  that  in  high-spin  state  m-LiMn02,  the  JT  distorted 
form,  is  more  stable  over  r-LiMn02,  which  agrees  well  with 
the  experiments.  While  for  low-spin  state,  the  changes  caused 
by  JT  distortion  are  almost  negligible  (see  Table  2),  which  con¬ 
sists  with  the  little  structure  differences  and  the  comparable  total 
energies  between  these  two  symmetries  and  thereby  supports  the 
above  conclusion  that  it’s  possible  to  synthesize  r-LiMnC>2  by 
compressing  the  Mn— O  bond  to  keep  the  Mn3+  ion  in  low-spin 
state. 

Based  on  the  aligned  total  DOS  (Fig.  3),  we  can  calculate 
the  change  with  JT  distortion  in  the  total  electronic  energy  per 
formula  unit  (A£eiec)  by  evaluating  Eq.  (3)  [6].  The  limits  of 
integration  are  from  the  lowest  Li- Is  state  to  the  E?.  Values 
of  —3.876  and  — 0.640  eV  are  calculated  for  high-spin  and 
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Table  2 


The  net  charges  of  Li,  Mn  and  0  ions,  bond  overlap  populations  (OP)  of  Li— O,  Mn— O  bonds  and  the  corresponding  bond  lengths  ( L )  before  (r-LiMn02)  and  after 
(m-LiMn02)  JT  distortion  for  high-spin  and  low-spin  states 


SflX 

Cell 

Li 

Mn 

O 

OPii-o 

Lu-o  (A) 

OPMn-0 

o 

^Mn-0  (A) 

4 

r-LiMn02 

+1.03 

+0.61 

-0.82 

-0.02 

2.1879 

+0.88 

2.0665 

m-LiMn02 

+1.04 

+0.58 

-0.81 

-0.04 

2.1008 

+0.90 

1.9679 

2 

r-LiMn02 

+1.05 

+0.42 

-0.73 

-0.06 

2.1086 

+0.98 

2.0021 

m-LiMn02 

+1.06 

+0.41 

-0.73 

-0.06 

2.0929 

+0.98 

2.0014 

low- spin  states,  respectively.  In  the  simplest  description  for 
the  JT  distortion,  the  JT  stabilization  energy  consists  of  two 
parts:  an  electronic  term  that  is  linear  in  the  distortion  (5) 
and  an  elastic  term  that  varies  quadratic  with  5,  which  can 
be  written  as  Eq.  (4)  [6,22],  where  A  is  the  electron-phonon 
coupling  parameter  and  k  is  an  elastic  constant.  Roughly,  Ejt 
equals  to  the  decrease  in  the  total  energies  A E.  Therefore,  the 
increase  in  the  elastic  energies,  i.e.  the  second  term  in  Eq.  (4) 
can  be  calculated  by  the  difference  of  A E  and  ALeiec-  With  the 
definition  of  8  as  the  difference  between  the  longest  and  shortest 
Mn— O  distances  in  the  elongated  octahedron,  we  can  obtain 
A  = -9.623  eV/A,  /c  =  45.622  eV/ A2  and  A.  = -711.11  leV/A, 

z  °  O 

k  =1.560  x  10  eV/A“  for  high- spin  m-LiMn02  and  low- spin 
m-LiMn02.  Obviously,  the  later  k  value  is  too  large  to  be  rea¬ 
sonable,  which  reinforces  the  conclusion  again  that  in  low-spin 
state,  the  JT  distortion  should  not  take  place  and  the  material 
should  have  the  R3m  symmetry.  In  another  words,  r-LiMn02 
should  be  stabilized  over  m-LiMn02  in  low-spin  state. 

Eeiec  =  J  S  X  N(e)  ds  (3) 

E]j  =  X8  +  - k8 2  (4) 

4.  Conclusions 

First-principles  calculations  in  the  GGS  have  been  employed 
to  investigate  the  JT  effect  in  layered  LiMn02,  r-LiMn02  and 
its  JT  distorted  form  ra-LiMn02,  with  different  spin  configu¬ 
rations.  It  is  found  that  non- spin  state  is  instable  for  both  the 
symmetries.  In  low-spin  state,  JT  effect  should  be  absent  and 
the  layered  structure  should  be  stabilized  in  r-LiMn02,  which 
may  be  realized  by  managing  to  compress  the  Mn— O  bond. 
Only  in  high-spin  state,  JT  effect  is  active  with  an  energy  ben¬ 
efit  of  175  me V  and  the  layered  structure  adopts  the  form  of 
m-LiMn02,  which  agrees  well  with  the  experiment.  Through  JT 
distortion,  regular  Mn06  octahedron  has  been  elongated  with  a 
degree  as  large  as  20.02%,  and  thereby  it  is  predicted  that  the  JT 
distortion  will  vanish  at  1910  K.  Simultaneously,  JT  distortion 
has  compressed  the  structure  perpendicularly  to  the  triangular 
layers  but  relaxed  in  the  layers  as  well  as  a  slight  expansion 
of  0.54%  in  the  volume,  which  directs  useful  ways  to  suppress 
the  JT  distortion.  JT  distortion  also  leads  to  distinct  changes  in 
the  electronic  structures,  where  the  initially  degenerate  eg  bands 
have  been  open  by  a  direct  gap  of  687  meV  and  the  ligand  field 
splitting  has  decreased  remarkably.  In  chemical  bondings,  JT 


effect  has  increased  the  ionicities  of  Li  ion  but  decreased  that  of 
Mn  and  O  ions.  The  ionic  interaction  in  Li— O  bonds  as  well  as 
the  covalent  interaction  in  Mn— O  bonds  has  been  enhanced,  and 
the  participation  of  Li  ion  in  the  JT  effect  is  detected.  At  last, 
the  JT  stabilization  energy  of  m-LiMnC>2  has  decomposed  into 
an  electronic  term  of  3.876  eV  and  an  elastic  term  of  —3.701  eV, 
where  the  corresponding  electron-phonon  coupling  parameter 

o 

and  elastic  constant  has  been  calculated  ab  initio  as  —9.623  eV/A 
and  45.622  eV/A  ,  respectively. 
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